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We present a relativistic chiral effective theory for symmetric and asymmetric nuclear matter 
taken in the Hartree-Fock scheme. The nuclear binding is insured by a background chiral invariant 
scalar field associated with the radial fluctuations of the chiral quark condensate. Nuclear matter 
saturation is obtained once the scalar response of the nucleon generating three-body repulsive forces 
is incorporated. For these parameters related to the scalar sector and quark confinement mechanism 
inside the nucleon we make use of an analysis of lattice results on the nucleon mass evolution with 
the quark mass. The other parameters are constrained as most as possible by standard hadron and 
00 ■ nuclear phenomenology. Special attention is paid to the treatment of the propagation of the scalar 

' fluctuations. The rearrangement terms associated with in-medium modified mass and coupling 

constants are explicitly included to satisfy the Hugenholtz -Van Hove theorem. We point out the 
, important role of the tensor piece of the rho exchange Fock term to reproduce the asymmetry energy 

5_( ■ of nuclear matter. We also discuss the isospin dependence of the Landau nucleon effective mass. 
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(N ! I. INTRODUCTION 

A fundamental question of present day nuclear physics is to relate low energy non perturbative QCD and 
"^T* ' in first rank chiral symmetry and confinement to the very rich structure of the nuclear many-body problem. 

. However it is presently hopeless to derive the observed nuclear properties from the underlying QCD and a 
^ ' more modest ambition is to put some constraints on the modelling of nuclear matter properties not only from 
'. hadronic phenomenology but also from lattice QCD data. This may also constitute a starting point to elucidate 
an old and central question of strong interaction physics, namely the interrelation between the many-body effects 
| governing the equation of state of nuclear matter and the nucleon substructure response to the nuclear environmment. 

Q\ ' A first attempt to go beyond the standard non relativistic treatment of nuclear matter is the relativistic mean field 
t-H , approach initiated by Walecka and collaborators [l[ . In this framework the nucleons move in an attractive background 
scalar field and in a repulsive vector background field. This provides a very economical saturation mechanism and 
a spectacular well known success is the correct magnitude of the spin-orbit potential since the large vector and 
C*~) scalar fields contribute to it in an additive way. Another successful modern attempt is based on in-medium chiral 
| perturbation theory where a pion loop expansion is performed on top of scalar and vector background scalar fields 
using a density functional formulation Now the question of the very nature of these background fields has to 
be elucidated or said differently it is highly desirable to clarify their relationship with the QCD condensates and in 
particular the chiral quark condensate. 

To address this question we take the point of view that the effective theory has to be formulated, as a starting point, 
in term of the fields associated with the fluctuations of the chiral quark condensate parametrized in a matrix form 
W = <t + if ■ if . The sigma and the pion, associated with the amplitude and phase fluctuations of this condensate are 
promoted to the rank of effective degrees of freedom. Their dynamics are governed by an effective potential, V(a, w), 
having a typical mexican hat shape associated with a broken (chiral) symmetry of the QCD vacuum. Explicit 
construction of such an effective theory for the description of nuclear matter can be performed for instance within the 
NJL model [|. 

As proposed in a previous paper Q an alternative and very convenient formulation of the resulting sigma model is 
obtained by going from cartesian to polar coordinates i.e., going from a linear to a non linear representation, according 
to: W = <T + if-n = SU = + s) exp (if ■ (p^jj^). The new pion field <p„ corresponds to an orthoradial soft 
mode which is automatically massless (in the absence of explicit chiral symmetry breaking) since it is associated with 
rotations on the chiral circle without cost of energy. The new sigma meson field, S, which is a chiral invariant, describes 
a radial mode associated with the fluctuations of the "chiral radius" around its vacuum expectation value, f^. It 
can be associated with the ordinary sigma meson which gets a very large width from its strong decay into two pions. 
Since it has derivative couplings to the pion field, it decouples from low energy pions whose dynamics is described 
by chiral perturbation theory. The evolution of the expectation value of S is related to the non pionic contribution 
to the in- medium chiral condensate [U, [1[ . This expectation value plays the role of a chiral order parameter around 
the minimum of the effective potential and the medium can be seen as a shifted vacuum. With increasing density, its 
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fluctuations s = S — f„ are associated with the progressive shrinking of the chiral circle and it governs the evolution 
of the nucleon mass. Here comes our main physical assumption proposed in ref. We identify this chiral invariant 
s field with the sigma meson of nuclear physics and relativistic theories of the Walecka type, or, said differently, with 
the background attractive scalar field at the origin of the nuclear binding. This also gives a plausible answer to the 
long-standing problem of the chiral status of Walecka theories. 

One motivation of the present work is to study in some details whether this hypothesis yields a viable description 
of nuclear matter. It is nevertheless well known that in such chiral theories, independently of the details of the 
modelling, tadpole diagrams associated with the mexican hat potential automatically generate attractive three-body 
forces destroying saturation 0, Q . The origin of this failure can be attributed to the neglect of the effect of nucleon 
substructure linked to the confinement mechanism as already pointed out in some of our previous works 0, 0, [H • 

Our article is organized as follows. The second section, which is in some sense a brief summary of our previous 
works, is devoted to the constraints brought by lattice data for the description of the nucleon and nuclear matter. 
In section IIIII we present the chiral lagrangian and section IIVI is devoted to the construction of the hamiltonian in 
the static approximation; we also give a detailed description of the treatment of the propagation of the in-medium 
modified scalar field. The Hartree-Fock approach including rearrangement terms is presented in section [V] and is 
applied to the case of infinite matter in section IVII Finally in section IVIII numerical results are given and the results 
discussed. 



II. CONSTRAINTS ON THE CHIRAL EFFECTIVE THEORY FROM QCD SUSCEPTIBILITIES 

A. Tests of the effective theory with a chiral invariant scalar field 

Once the appropriate couplings of the chiral fields to the baryons are introduced one can build an effective lagrangian 
to describe nuclear matter. Vector mesons (a; and p) must be also included to get the needed short range repulsion 
and asymmetry properties (see sections IIIII and IVII[) . At the Hartree level, the pion and the rho do not contribute for 
symmetric nuclear matter whose energy density written as a function of the order parameter s = (s) is : 

T = £ ° = / 6{PF p) E ^ + V{5) + p2 - (1) 

e;{s) = v^TlfM is the energy of an effective nucleon with the effective Dirac mass M N (s) = Mm + gs s - gs is 
the scalar coupling constant of the model; in the pure linear sigma model it is gs = Mpf/fw The effective potential 
V(cr, 7?) when reexpressed in term of the new polar representation has the typical form : 

vM-in* (« 2 + 5j +•••)• 

s is obtained by minimization of the energy density and is given at low density by : s ~ —(gs/rn^) ps- Its negative 
value is at the origin of the binding but the presence of the s 3 term (tadpole) has very important consequences as 
already mentioned in the introduction. This tadpole is at the origin of the chiral dropping [9[ of the sigma mass 
Am* ~ — (33s/2/tt)ps (a ~ 30% effect at po) and generates an attractive three-body force which makes nuclear 
matter collapse and destroys the Walecka saturation mechanism. Hence the chiral theory does not pass the nuclear 
matter stability test. 

This failure, which is in fact a long-standing problem 0, 0], is maybe not so surprising since the theory, as it is, 
also fails to describe some nucleon structure aspects as discussed below. The nucleon mass, as well as other intrinsic 
properties of the nucleon (sigma term, chiral susceptibilities), are QCD quantities which are in principle obtainable 
from lattice simulations. The problem is that lattice calculations of this kind are not feasible for quark masses smaller 
than 50 MeV, or equivalently pion mass smaller than 400 MeV, using the GOR relation. Hence one needs a technics 
to extrapolate the lattice data to the physical region. The difficulty of the extrapolation is linked to the non analytical 
behaviour of the nucleon mass as a function of m q (or equivalently m~.) which comes from the pion cloud contribution. 
The idea of Thomas et al [lfl ] was to separate the pion cloud self-energy, E n (m n , A), from the rest of the nucleon mass 
and to calculate it in a chiral model with one adjustable cutoff parameter A. They expanded the remaining part in 
terms of as follows : 

M N (ml) = a + a 2 ml + a±m% + E-^m^A) . (2) 

At this point it is important to stress that the above expansion is in reality an expansion in terms of the current 
quark mass m q which is the genuine parameter occuring in the lattice calculation. The pion mass appearing in 
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eq. <j2j) is just the pion mass deduced from the quark mass assuming the GOR relation m 2 = —2m q (qq) vac /f 2 . 
This reparamctrization has been adopted only for convenience. The best fit value of the parameter 04 shows little 
sensitivity to the shape of the form factor, with a value 04 ~ — 0.5GcV~ 3 while 02 ~ 1.5 GeV -1 (see ref. The 
small value of 04 reflects the fact that the non pionic contribution to the nucleon mass is almost linear in (i.e., in 
m q ). Taking successive derivatives of Mn with respect to (i.e., to m q ), it is possible to obtain some fundamental 
chiral properties of the nucleon, namely the pion-nucleon sigma term and the scalar susceptibility of the nucleon. The 
non pionic pieces of these quantities are given by : 



non— pion UJ.YJ.pj 2 uivj 2 i o 4 oni\iT \t I o\ 

t n =m q g— - m *g^2 =fl 2< + 2ffl 4<- 29 MeV - ( 3 ) 



o / non — pion /o \ / \ 

Xns ~ ^ - ^Tdmi\-^—)-^r A - () 



In the above equations the first equalities correspond to the definitions, the second equalities make use of the GOR 
relation and the last ones come from the lattice QCD analysis. With typical cutoff used in this analysis, A ~ 1 GeV, 
which yields a N ~ 20 MeV, the total value of the sigma term is crpj ~ 50 MeV, a quite satisfactory result in view 
of the most recent analysis. It is interesting to compare what comes out from the lattice approach with our chiral 
effective model. At this stage the only non pionic contribution to the nucleon mass comes from the scalar field, or 
more microscopically the nucleon mass entirely comes from the chiral condensate since the nucleon is just made of 
three constituent quarks with mass Mq = g(S) vac = gf„ ~ 350 MeV. Hence the results for the non pionic sigma term 
and scalar susceptibility arc identical to those of the linear sigma model : 

<r N -I*9s m 2, X NS - * % m 4 ■ ( b > 

The identification of a ™ n ~P lon w ith crffi of our model fixes the sigma mass to a value m a = 800 MeV, close to the one 
~ 750 MeV that we have used in a previous article ■ As it is the ratio gs / m 2 which is thus determined this value of 
m a is associated with the coupling constant of the linear sigma model gs = Mpj/f^ = 10. Lowering gg reduces rn a . 
Similarly, the identification of Xns wl ^ [i the lattice expression provides a model value for 04. The numerical result is 
= —3.4 GeV" 3 while the value obtained in the expansion is only —0.5 GeV -3 . 



B. Nucleon structure effects and confinement mechanism 



The net conclusion of the above discussion is that the model as such fails to pass the QCD test since the 04 coefficient 
is much larger in the chiral model than the one extracted from the lattice data analysis. In fact this is to be expected 
and even gratifying because it also fails the nuclear physics test as discussed in subsection III A[ We will see that 
these two important failures may have a common origin. Indeed an important effect is missing, namely the scalar 
response of the nucleon, kns = d 2 Mpi /ds , to the scalar nuclear field, which is the basis of the quark-meson coupling 
model (QMC) introduced in ref. [ll|. The physical reason is very easy to understand: the nuclcons are quite large 
composite systems of quarks and gluons and they should respond to the nuclear environment, i.e., to the background 
nuclear scalar fields. This response originates from the quark wave function modification in the nuclear field and will 
obviously depend on the confinement mechanism. This confinement effect is expected to generate a positive scalar 
response kns, i- e -, it opposes an increase of the scalar field, a feature confirmed by the lattice analysis (see below). 
This polarization of the nucleon is accounted for by the phenomenological introduction of the scalar nucleon response, 
knsi hi the nucleon mass evolution as follows : 

M N (s) = M N + g S s + 1;k N ss 2 + .... (6) 

This constitutes the only change in the expression of the energy density (eq. [I} but this has numerous consequences. 
The effective scalar coupling constant drops with increasing density but the sigma mass gets stabilized : 

dM% M N _ * 2 d 2 e 2 /3 5s v 
3s( s ) = —g=- = -J 1~ k ns s i m a = — ~m a -(— Kns)Ps- (7) 
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The non-pionic contribution to the nuclcon susceptibility is modified, as well [8| : 

(a) (qq)lac ( 1 1 \ 1 o {m)% ac 1 (3gs \ , a , 
Xns = {^-KJp = K {IT - KNS ) • (8) 

We see that the effect of confinement (kns) is to compensate the pure scalar term. Again comparing with the lattice 
expression one gets a model value for the 04 parameter : 

<U = -^(3-2C). (9) 

where C is the dimensionless parameter C =(J^J2 M) kns- Numerically 04 = —0.5 GeV~ s gives C = +1.25, implying 
a large cancellation. As discussed in ref. 0, @, Q such a significant scalar response will generate other repulsive forces 
which restore the saturation mechanism. At this point it is important to come again to the underlying physical picture 
implying that the nucleon mass originates both from the coupling to condensate and from confinement. In the original 
formulation of the quark coupling model, nuclear matter is represented as a collection of (MIT) bags seen as bubbles 
of pcrturbative vacuum in which quarks are confined. Thus in such a picture the mesons should not appear inside the 
bag and should not couple to quarks as in the true non perturbative QCD vacuum. Consequently the bag picture is at 
best an effective realisation of confinement which must not to be taken too literally. Indeed, QCD lattice simulations 
strongly suggest that a more realistic picture is closer to a Y shaped color string (confinement aspect) attached to 
quarks [l2T ]. Outside this relatively thin string one has the ordinary non perturbative QCD vacuum possessing a chiral 
condensate from which the quarks get their constituent mass. 



III. THE CHIRAL LAGRANGIAN 



To get a more complete description of symmetric and asymmetric nuclear matter we complete the lagrangian used 
in ref. Q essentially by the introduction of the rho meson. We also consider the possibility of incorporating the scalar 
isovector 5 meson. Written with obvious notations, it has the form 



with 



C = #*t"3 m ¥ + C s + £ w + C p + C s + C v (10) 

C s = -M N (s)^ - V(s) + ±d»8d M s 

Cu> = -&,w„*7"tf + ^u"w„ - \ f» v F^ 

C P = -g p p a p^T a ^ - g p J^-d vPa ^a^T a ^ + - m 2 p pa^fc - ^G^G a p„ 

Cs = -g s 6 a *T a * -\m 2 s 5 2 + \d»5d^5 

Ar = ^-3^fT"7 5 r a *- \ ™l<fiL + ^'Va^¥W (11) 

As discussed previously the form of M^v(s) (eq. [5]) reflects the internal nucleon structure and contains a quadratic 
term involving the scalar response of the nuclcon which is constrained by lattice data. However the nucleon mass may 
very well have higher order derivatives with respect to the scalar field. In practice, as in our previous works 0,13 1 we 
introduce a cubic term : 

M N (s) = M N + gsS + I KNS (s 2 + ^j-^ . (12) 

Hence the scalar susceptibility becomes density dependent 

8 2 M N ( s\ 

and vanishes at full restoration, s = —fn, where s is the expectation value of the s field. Hidden in the above 
Lagrangian is the explicit chiral symmetry breaking piece 

C x sb = co = --Tr(U + s)exp{if-(p n /f v ) ~ cs - Y^fl ( 14 ) 
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which generates the pion mass term with the identification c = f„ rri^. It is thus implicit that neglecting the higher 
order terms in the exponent, the self-interactions of the pions are omitted. Notice that the only meson having a 
self-interacting potential V(s) is the scalar meson s. We take it in practice as in the linear sigma model with the 
inclusion of the explicit chiral symmetry breaking piece : 



V(s) 



A 



((/* + sf -J) 2 ~ Umls 



2/* 



8/ 2 



(15) 



The other parameters {g LJl g p , n p , qa and the meson masses) will be fixed as most as possible by hadron phenomenology 
It is in principle also possible to calculate or at least to constrain these parameters in an underlying NJL model. One 
specific comment is in order for the tensor coupling of vector mesons. The pure Vector Dominance picture (VDM) 
implies the identification of n p with the anomalous part of the isovector magnetic moment of the nucleon, i.e., n p = 3.7. 
However pion-nucleon scattering data [l3| suggest k p = 6.6 (strong rho scenario). We will come to this point later on 
in the discussion of the results ( section IVII[) . The omega meson should also possess a tensor coupling but, according 
to VDM the corresponding anomalous isoscalar magnetic moment is = 0.13. Since it is very small we neglect it 
here. For completeness we also add a delta meson which may generate a splitting between the proton and neutron 
masses but with the chosen coupling constant, gs — 1, its influence is in practice negligible. 



IV. CONSTRUCTION OF THE HAMILTONIAN 



The conjuguate momenta of the various mesonic fields are 



n 



F 0j , 



n„ 



d Sa 



H ajp = G^ + g p ^JL-Hfo> U T a V, 



n 



A Jtt 



(16) 



The hamiltonian is obtained by the usual generalized Legcndre transformation with the result : 



H = dr 



* -if . V * + H s + H u + H p + H s + H v 



(17) 



with : 



H 5 



dr 
dr 
dr 

2 [ 9 > 



M N (s)W + | (n 2 . + (Vs) 2 ) + V(s) 



g^ * 7 ^* + i ( (jivf - m 2 uj^ujp + vw 3 ■ w - (v • ay 



gpPan^Ta^f + g p7 T^- I d jPai ^a ij T a ^ - ll al *cr ?0 T Q * 



£ I {ll ap f - mlp» Pap + Vpi ■ Wpi - (V • p a f I - Ii op • V Pa 



gsSa^ + - (n 2 5 + (V<5 Q ) 2 + ml8 2 a 



9A 

2/. 



*7 5 7°r a *) 



(18) 
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A. Static approximation 



We will first formulate the Hartree-Fock approach in the static approximation, i.e., neglecting the retardation 
effects. In such a case the conjuguate momenta are : 

n s = o, it = VuA u aS = o, u ajp = v jP ° a + 9 P ^tt *^°ra*, n oj - * 7 5 7°ra*. 

I M N I J n 

The various pieces of the hamiltonian simplify according to : 



jjstatic 
jjstatic 



jjstatic 



dr 

dr 

f 
2 

= / dr 



M N (s)M + i(Vs) 2 + V(s) 
9u> 0Jn *7 M * - \ (ml ui^ + Vw" • Vw^ + (V • u)' 

l IVljq 



^nP^Pau + V^-Vp aM + (V-p a ) 



Hi 



This hamiltonian can be rewritten as 



= / dr 



gs 5 a *r a * + - ( (V5 a y + m£<£ 



dr ( K 



mesons ) 



The first term is the nucleonic piece including the Yukawa coupling of the nucleons to the meson fields 
K = V + M N (s) + g^oJul^ + g P p a u^T a + 9pY^d.jp ail a M T a 

„ gA ^ _ \ 

+ gS a T a + -— VfaTT ■ 7 l T a I *. 

The second piece os of purely mesonic nature : 

1 



H 



mesons 



2 (VS)^ + V(S) - | ( m 2 ^ + VcjM . V ^ + (y . 



1 



+ l((VS a y + m ffi) + \ ((V^w) 2 + mlvl 



(19) 
(20) 



(21) 



(22) 



B. Equation of motion for classical and fluctuating meson fields 



In the static approximation the equation of motion for each meson field ipn can be written formally as [H, II#] = 0, 
where Yin is the conjugate momentum of the field ipa. This gives : 

dK dM N 
ds 



-V 2 s + V'(s) 



dK 



-V 2 o^ + m 2 + ^ift(V ■ (3) = = *7 M * 
~V 2 p^ + m 2 p» + <^(V ■ Pa) = ^- = 9p ^r a * g p ^- d 3 (*o«V *) 

V^ a + 777| $ = T^T - = g8 *T a * 



2 i 2 

'" a7r + 777 w V? a7r 



-VV 



d^ajr 2 f n 



V • ^7 5 7T a *. 



(23) 
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Following the method used in ref. [14| we now assume that it makes senses to decompose each meson field as : 

Vr = <Pr. + A(p R (24) 

where (p^ ; = (</?r) denotes the ground state expectation value of the meson field tpa and Aip^ corresponds to its 
fluctuation considered as a small quantity. 

Since the case of the scalar field is the most delicate one we treat it below with some details. The equation of 
motion for the s field can be expanded in As according to : 

BK 8 2 K 

- V 2 (s + As) + V'(s) + As V"(s) = -—(s) - As -5-5- (a). (25) 

us os z 

Explicitly we have : 

dK, , dK . - _ . dM N {s) 

— {s) = —=g%W with g* s = g _} =g s + n NS s + ... 

3 2 i^ . d 2 if . - . , . 9 2 A/ W (s 2 ) 

— - s = — - = kns with k WS = 7^-^ = kns + •■■ (26) 

os £ os z as 

We now develop the source term of the equation of motion according to : 



as \ os I \ os J \ os I \ as \ as / 



(27) 



and we consider the fluctuating term A(dK/ds) as small and of same order than As. In the same spirit we replace 
the second derivative of K by its expectation value : 

The equation of motion will be solved order by order : 

-VH + V'(s) = -(^=- g * s <**) 

-V 2 (As) + mf As = - (2f - (H)) = -9s (*« - <**» . (29) 

In the equation for the fluctuating part it appears the effective scalar mass : 

m* a 2 = V"{s) + k NS (30) 

already introduced in our previous work Q. This is the physical in-medium scalar mass propagating the quantum 
fluctuations of the scalar field, i.e., the quantum fluctuations of the chiral condensate in a non-trivial way. In 
particular we will see below that it is the mass appearing in the Fock term of the scalar exchange at variance with 
V'(s)/s appearing in the Hartree scalar exchange. In that sense the treatment of the self- interacting scalar field 
deviates from the one of ref. UM ■ 



C. Kinetic, Hartree and exchange hamiltonians 



We now develop the hamiltonian (eq. I20[) to second order in the fluctuations, limiting ourselves to the nucleon field 
and the s field, the generalization to the other mesons being straightforward. 



H s = dr 



dK 1 , . n2 d 2 K 1 , 

1 - I — sV s - 2 AsV 2 s 



K(s) + As— + -JAs) ^ 



+ V(s) + AsV'(s) + -(As) 2 l/"(s) 



(V(As)) ; 



(31) 



Using the classical equation of motion and the one for the fluctuating field and replacing again the second derivative 
of K by its expectation value, we obtain : 



Hs = dr 



-Z7- V + M N (s)) * + 1 (V(s)) 2 + V(s) + l-g* s (** - As 



(32) 
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The approach can be extended to the other mesons. Only the delta field and the time component of the rho and 
omega mesons have a non zero expectation value, solution of the classical equations : 



-V 2 w° + ml w° = g u 



N 



-V 2 <5 a + m 2 a S = gs (*r a *) . 
The fluctuating fields are solutions of : 

-V 2 (A^) + ml Aw" = P^g^ (#7"* - <*7"*» 
-V 2 (Ap£) + m 2 Ap" = P»g p (*7*V tf - (#7"T a tt) 



2M N 



8 3 (Va vj T a V) 



2M N 



-V 2 (A<5 a ) 



A<5 a =g s (*r a * - (*r a *» 



-V 2 (A<^ aT ) + to 2 A<p Q 



2/. 



with 



p o 



1, P° = 0, P l = 0, /'; ; /",(.,•; 



(33) 



(34) 



(35) 



Generalizing the result of eq. (|32p to all the mesons, the full hamiltonian, in the static approximation, can be written 
in a form reminiscent of the density functional theory : 



H H-kin+Hartree ~t~ H xc 



(36) 



The first term is a one-body operator containing the kinetic energy hamiltonian of the nucleons and the other pieces 
of the hamiltonian contributing to the Hartree energy. Its explicit form is : 



H, 



kin-j- Hartree 



dr 



* - ij ■ V 



M N (s) + g w ^°7o + g P plloT3 



■9p 



2M N 



Vis) + \{Vs) 2 -\mliwy -\(y^) 2 



1 



■0\2 



m (Pi) 



1 



m 2 & Si + 



1 



V<5 ; 



(37) 



Considering only this part of the hamiltonian we come to the conclusion that symmetric and asymmetric nuclear 
matter are seen as an assembly of nucleons, i.e., of Y shaped color strings with massive constituent quarks at the 
end getting their mass from the chiral condensate. This nucleons move in self-consistent scalar (s,Ss) and vector 
background fields (o) ,p 3 ). The scalar field which we associate with the radial mode of the condensate modifies the 
nucleon mass according to Mjv(s) = Mjv + gs s + knss 2 /2 + .. . gs s describes the in-medium dropping of the 
nucleon mass from its coupling to the in- medium modified chiral condensate but there is another term, kjvs s 2 , which 
corresponds to the response of the nucleon to the background scalar field. The scalar response kns depends on the 
structure of the nucleon and takes into account the modification of the quark wave functions inside the nucleon and 
obviously depends on the confinement mechanism. The nucleon is also submitted to an isovector field, (Q ) and in 
asymmetric nuclear matter to an isovector vector field (p®) or even an isovector scalar field (S3). 

The second piece of the hamiltonian, H xc incorporates the exchange term mediated by the propagation of the 
fluctuations of the meson fields. In particular the scalar fluctuation, i.e., the fluctuation of the chiral condensate 
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propagates, as already stated, with an in-medium modified sigma mass to* 2 = V"(s) + k-ns Ps- Its explicit form is : 

f 



H ,■ 



dr ■ 



g* s As A + g u Au, L (tf 7 "tf ) 



- 9 P 



gs AS a A (*r a *) 



9A 

2fn 



where we have used systematically the notation : 

A(*r*) = *r<£ - <#r*) . 

We introduce the (static) propagators for the fluctuating fields : 

(-V 2 + m?(r)) D a {r - r') = <5< 3) (r - r') 

(-V 2 + to 2 ) D a (r - r') =<5 (3) (r - r'), a =w,p,S,% 

and solve formally for the fluctuating fields. H xc can be written as : 
H xc = - I drdr' 



- g* s (r) ffs(r') A (**) (r) D a (r - r') A (M) (r') 

gl A (*y*¥) ( r ) D ufiv (v - r') A (§7"*) (r') 
5 2 A (*yV a tf) (r) D p ^(r - r') A (*7 tf r a *) (r') 
2 5 ?^- A (r) a,-D p ^(r - r') A (*t"t *) (r') 



2Mjv 



A (*<r"V 4-) (r)a,5p p ^(r - r') A (*^'r„*) (r') 



- 5 2 A (*r Q *) (r) D 4 (r - r') A (*r a *) (r') 
\ 2 

Yf ) (*7 6 7V,#) (r) fi^A^r - r') (fryVVa*) (r') 



(38) 



(39) 



(40) 



where we have introduced the tensor propagator D^^ v (v — r') whose non vanishing components are : D^oo(r — r') = 
£> w (r — r') and D^ijir — r') = (Sij — didj/m^) D^i? — r') and a similar one for the rho meson. 



V. HARTREE-FOCK APPROACH 



A. Hartree-Fock energy 



In the Hartree-Fock approximation the ground state is represented by a Slater determinant made of single particle 
states with Dirac wave functions <p^ {r)XN, with N = p,n for protons and neutrons. The various densities which 
appear in the sources of the classical equations of motions are : 



E ^< + = PSp + PSn 

E- V p — n n 

Va^Pa ~ { Pa { Pa= PSp - PSn 

E - ^>°v^ 



(41) 



a<F 
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The kinetic plus Hartrcc piece of the total energy is : 

Ekin+H 



J dr J2 ^H7 • V + M N (s) + g u w° 7o + g P f%Jo 



a<F 



'9p 



2 AT 



N 



Vp°-I] 75 + g s S 3 



[dr\Y, ^H7"V + M N (s) + g 

** ^ TP* 



• u uj° 70 - g p pi 70 



a<F 



+ Jdr + \ (v.") 2 - ^m 2 , (^) 2 - i (Va>° 

- \ (Vp°) 2 + im^l + \ (W 3 ) 2 
The Fock term contribution to the energy comes entirely from H xc . It can be split as : 



i'Focfe - ^jr ock + i^Fock + ^Fock 



P (<5) , 

^Fock ^Fock 



We give here the expression for the scalar field piece 



E 



0) 



_Focfe 



= ~ j dvdv' Tr(S p (v' - r) S p (r - r') + S„(r' - r) S n (r - r')) 
^(r)<£(r')Ar(r-r'), 

the other contribution being given in appendix A. Here S p and S n are matrices in Dirac space : 

(S p (v-r')) a(j = £ (*>S) a (r) (^)^(r') (^n(r-r'))^ = £ (r) (r'). 



(42) 



(43) 



(44) 



a<F 



a<F 



B. Hartree-Fock equations 

The single particle orbitals are obtained by minimizing the HF energy with respect to the (pP' n (x), with the 
constraint that the single particle wave functions are normalized, i.e., 



= 



(45) 



and the Lagrange parameters have to be identified with the single-particle energies. According to the results of 
the previous section the structure of this HF energy is : 



E = 



a,N J 



Va ( - i7 ■ V + M N (s) + VrTr + i h ( P V<p R -e R )^ 



(NN 

9rr 



Y ^f— drdr'Tr(S N (r'-r)T R S Nr (v-r')T R ,)D RRI (r-v') 

3/ . AT AT/ J / 



(46) 



R,R';N,N 



The first two lines correspond to the kinetic plus Hartree energy (eq. I42p . the second line being the purely mesonic 
piece; the last line represents the Fock (exchange) term (eq. S3])- The various coupling constants, g^ R \ , \ 
Dirac operators, T R ,& R , and propagators, D RR >(r — r'), are obtained in an obvious way by direct comparison with 
the explicit expressions of the energy given above (eqs. |43l l44l and Appendix A). For the the particular case 

of the scalar channel the coupling g^J? has to be inserted inside the integral and replaced by <7<j(r) gg(r'). In the 
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minimization procedure the derivative of the expectation values of the mesonic fields tp R is not taken since it is already 
accounted for by the equations of motion. The last line of the HF energy actually generates two kinds of contributions 
in the HF equations; the first one comes from the functional derivative of the nucleon propagators, 5jv(r — r') and the 
second one comes from the derivative of g$ g$ D a in the scalar channel (eq. I44[) . The later is actually a rearrangement 
term associated with the many body forces originating from the in-medium dependence of the scalar meson mass 
and of the scalar coupling constant. Ignoring for the moment the rearrangement terms, the Hartree-Fock equations 
writes : 



dx' (r 



h N 

a (ord) 



r'} ^(r')=«(r) 



(47) 



where the single particle hamiltonian, h N , which is represented by a matrix in Dirac space, has a local piece (the 
direct Hartree piece) and a non local piece (the Fock exchange term). With the previous schematic notations it is 
such that : 



uN 
(ord) 



V + M N (s) + Y^g^cpR^Tn + ih^Vmir) ■ Q R )s^(v - r') 



R,R':N> 



(NN') 



(NN' 

9rr- 



E ^-D RR ,{r'-r) r ff S w -(r-r')r 



R,R':N> 



(48) 



Let us now come to the rearrangement term originating from the density dependence of the scalar mass and coupling 
constant. Its contribution to the single particle hamiltonian is defined by : 



uN 

(rg) 



fa) = \J rfxdx' ^T7-(Sm(x - X')S M (X' - X ; 



mi?) 



g* s (x)g* s {x')D a {x - x') 



To get the functional derivative of the scalar coupling constant we start from 

- kns (x) 



and take the functional derivative of the equation determining s (first equation of eq. (|29p ) 

5s(x) 



-V£ + V"(s(x)) + k NS (x) Ps (x) 



8<p»(T 



(49) 



(50) 



(51) 



Notice that the full in-medium sigma mass and the full inverse sigma propagator D a 1 appear in the left-hand side of 
the above equation. The solution of this equation is thus : 



<Sa(x) 



= - / dx'A^x - x')^(x')^(r)J(x' - r) = - D a (x r) 3 *(r)^(r 



It follows that the derivative of the scalar coupling constant has the explicit form 



= -K NS (x)D a (x- r) 5 J(r)^ v (r). 



(52) 



(53) 



Similarly the functional derivative of the sigma propagator can be obtained by taking the functional derivative of eq. 
(l39l. One obtains : 



< 2 (x) 



6D a {x - x') <5m* (x) 



8<p»(T 



D a {x - x'). 



(54) 
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(c) 




(b) 




gm 




g* s (x) 



K N s(r) 



(e) 



FIG. 1: Rearrangement terms; the crosses indicate the opening of a line to generate the various contributions to the single- 
particle hamiltonian. Detailed explanations are given in the text. 



The functional derivative of m* (x)) = V"(s(x)) + kjvs(x) ps(x) is given by : 



6<p»(T 



+ k NS ( X )^(r) ^(x-r) 



+ (x) (- D a (x t) 9 Ut)<P?(t) 



+ Kjvs(x)^(r) <$( 3 >(x - r). 
The explicit form of the functional derivative of the in-medium sigma propagator follows 
8D a {* - x') f JT ,, Jrn* a (y) 



S<P?(r) 
dyD a (x - y) [V"'(s 



dyD a (x - y)"-^p±LD a {y - x') 
, Okns 



s<p2 ft- 
ps) (y)D a (y - x'Wy - r)<£(r) V *(r) 



£> ff (x - r)D a (r - x') k NS r) ^ (r). 



(55) 



(56) 
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The rearrangement single particle hamiltonian (cq |49[) can be decomposed in two terms. The first one comes from 
the functional derivative of <7g(x) : 

(r| 7 °/$ sl) |^) = -I J dxrfx'Trj/^x - x')5 M (x' - x)) £> ff (x - x') 

^ s (x)^(x')^(x-r)^(r)^(r) 
+ (k W s(x') 5 S(x)^(x' - r) 5 |(r)^(r) 
The second term originates from the derivative of the sigma propagator : 



(57) 



uN 
(rg2) 



<p%) = J dxdx , Tr M (s M (x - x')S M (x' - x)) <£(x) g* s (x' 



dyD c (x - y) [V"'(s 



dk NS 
ds 



Ps)(y)D a (y - x')A,(y - r) g* s (r) 



-£> CT (x - r)D a (r - x.')k ns t) 



Notice that these rearrangement terms are of local nature, hence the Hartree-Fock equations finally write 



dr' (r 



h N 



r') ^(r') + (r| 7 °^ 9l) |^) + 



uN 

rg2) 



7° h 



(58) 



(59) 



The diagrammatic interpretation of the various terms is straightforward as shown on fig. 1. Each contribution 
corresponds to the opening of the fcrmion lines in the Fock energy diagram. The ordinary non local piece of the 
single particle hamiltonian simply comes from the opening of the fermion lines (fig. la) ignoring the dressing of the 
sigma propagator and the in-medium modification of the scalar coupling constant. These medium effects appear in 
the other diagrams and are depicted schematically on figure lb,c,d,e. Opening the fermion loop renormalizing the 
scalar coupling constant (fig. lb) generates the first rearrangement term. Opening the fermion lines associated with 
the dressing of the sigma propagator (fig. lc,d,c) generates the second rearrangement term. 

The above formalism can in principle be directly used for finite nucleus calculations. However important specific 
difficulties occur which are due in particular to the density dependence of the sigma mass and to the non trivial 
structure of the rearrangement terms. In practice this requires more study and we limit ourselves in the present paper 
to infinite nuclear matter. 

VI. INFINITE MATTER 

In infinite nuclear matter, the single particle orbits arc plane waves, labeled by momentum and spin indices k = (k, s) 
for each isospin state N : 



«(k, s) e lk T XN 



The Fourier transform of the single particle potential can be introduced through 



(r^o^jr') = 



(2tt) 2 



s lk - (r - r ') 7 ° h N (k) 



(60) 



(61) 



(62) 



where Eg (k), Sy (k) and So (k) are the various self-energies of scalar or vector nature. Notice that the above 
quantities are a priori different for protons and neutrons. The Hartree-Fock equation becomes : 



and has always the general form 

7 ° h N (k) = M N + S s (k) + 7 • k( 1 + S y (k) )+ T »S„ (k) 



/i iV (k)u(k,s) =£>(k,s) 



(63) 
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This equation is formally identical to the Dirac equation for a free particle, giving for our normalisation (|60[) : 



2E* \ E * +M * X 

where 



u(k, g ) = A / -*tf* ( (X-k* „ ) 



£* = - S (k) = VM* 2 + k* 2 

Af* = Miv + S s (k) (65) 
k* = k (1 + £v (k)) 

are the effective energy, Dirac effective mass and effective momentum. Again notice that M* and E* are different 
for protons and neutrons. In the following formula we will distinguish them with p and n indices. We also introduce 
the Fermi momenta, pp for the protons and np for the neutrons, the occupation numbers for the protons, N p k — 
Q(pf — |k)|, and for the neutrons N n k = &(np — |k)| and the vector and scalar density for protons and neutrons : 

f 2dk [ 2dk 

Pp = /f, x 3 "pki Pn= 77. x 3 Nnk, p = p p + p n 



(2tt) 3 r " J (2tt) 3 

2dk Ar Ml [ 2dk Ar M* 

Ps P = I N p* ~e*> PSn = J J^y Nnk ~E*' ps = PSp + PSn - (66) 



The S p and S n matrices take the explicit form 



V ^ \ 2E 'off* 



Vpk 



(S n (r-r>)) a0 = I^ e *-(—') AU; (67) 

The energy density can be written as : 

^ ^kin+Hartree ^Fock 

(68) 

Once the equations of motion for the expectation value of the classical u>q, po and S fields have been used, the 
Kinctic+Hartrce contribution to the energy density has the form : 



£kin+ H artr ee 



2rfk / k* M*\ 

W?( l '? tl '" W F)/ VW 



The Fock contribution from scalar exchange is given by 



c 



(a) gf f dk dk' i v / M* M * k* k*\ 

Focfc 2 ,/ (2tt) 3 (2tt) 3 (k - k') 2 + mf ^ \ E* E'* E* E'* J k k ' V ; 



and the other Fock terms are listed in appendix B. 

Let us now come to the various contributions to the self-energy. Ignoring for the moment the rearrangement terms, 
they can be obtained starting from cq. (|48[) . For the protons the scalar component of the self-energy is : 

S s (k) = M N (s) - M N - % { PS p - PSn ) 



_ g<£P f dk' 1 1^ ( v r + M' - ]n 

+ 1 1n —J (k - k') 2 + m% A Tr [ Tr -2E*- Vr 'J N * 

g%P f dk' l ( v r + M'* r \ 



\r>* 
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the time-component part of the vector self-energy is : 

S (k) = 4" iPp + Pn) + ^(Pp~ Pn) 



\- g& f dk' 1 1^ ( r + \ „ 



2 J (2tt) 3 (k - k') 2 + m\ 4 V 2£;,: 
£ — / ( 27r) 3 (k _ kT + w2fl 4 J> (jo F fl ) AT k , (72) 



and the spatial part is 



9^' 

R,R';N 



dk' 1 1_ / - F + M 



(2tt) 3 (k - k') 2 + m| 4 V 



For the neutrons, the self-energies are the same as those for the protons except that the p and n indices are exchanged. 
Note that in infinite nuclear matter, the expectation value s does not depend on the position, therefore the Fourier 
transform of the s propagator is as simple as for the other mesons, simply replacing the vacuum mass by the in- 
medium sigma mass. Here we have decomposed the p term into a vector one, a tensor one and a crossed one between 
the vector and the tensor interactions. 

Finally, the rearrangement terms give momentum independent contributions to the scalar self energy : 

S S (rg i) - ~ 2 — e Fock (74) 



/ dk dk' f 1 

dps 2 J (2tt) 3 (2tt) 3 V (k - k') 2 + mf 



E* E* E* E 
n \ / N 

where the derivative of the in-medium sigma mass with respect to the scalar density is 

V (s) H — - p s - h k ns 



dp s ds , ) dps 

% (W) + - B 



*NS-4,(V'"(S) + ^PS). 



VII. RESULTS 



In the previous sections the results for the energy (eq. l70|) and the self-energy (eq. 171)72173]) have been presented 
in a way which exhibits the structure of the interactions and in particular the explicit form of the meson propagators. 
However in practice the angular integrations relative to the angle between k and k' are performed analytically and 
we obtain results which are similar in shape to the one used in ref. [l6| for the numerical calculations. 

As explained before the spirit of this paper is to study nuclear matter properties with parameters fixed as most 
as possible by lattice data and hadronic phenomenology. For the well established masses we take Mjy = 938.9 MeV, 
m u = 783 MeV, m s = 984.7 MeV, m T = 139.6 MeV and m p = 779 MeV. The pseudo-vector coupling constant of the 
pion is <m/(2/7t) where gA = 1-25 is the axial coupling constant and f n = 94 MeV is the pion decay constant. For the 
rho meson coupling constant we take the VDM value g p = 2.65. The scalar coupling constant in the vacuum is fixed to 
the linear sigma model value gs = 10. From a comparison with lattice data we deduced a sigma mass m a = 800 MeV. 
Hence the only parameters are the scalar suceptibility K^s (or the dimensionless parameter C = (/ 2 /2Mn)kns) and 
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the omega- nucleon coupling constant, g^. Nevertheless we do not take them as really free parameters but only allow 
small deviations around the lattice estimate C ~ 1.25 and the VDM/quark model value g u = 3 x 2.65 ~ 8. 

One uncertainty which is particularly important for the asymmetry energy, 05, is the tensor coupling constant of 
the rho to the nucleon. We will consider different cases between the pure VDM value k p = 3.7 and the value deduced 
from scattering data n p = 6.6 [l3l |. For each case the susceptibility parameter C and g u will be fixed. 

Finally we remark, as the Orsay group [HI, [l6[ , that the pion exchange and the tensor part of the rho exchange NN 
interaction contain a momentum independent piece which corresponds in configuration space to a contact interaction. 
In realistic many body calculations this S(r) contribution has to be suppressed by short-range correlations. The 
prescription of the Orsay group, also used in [Til , Ht| . was simply to remove the contact term, i.e., making the 
replacement in the Fock exchange terms : 



2 2 
q 2 q 2 



- 1 = - ■ 



q 2 + to 2 q 2 + m 2 q 2 + m 2 

2 ~2 ^2 



q q 



q 2 + m 2 q 2 + to 2 q 2 + m 2 ' 



m„ 

-i = -^-r^2- ( 77 ) 



In other words, the pion and tensor rho exchange become ordinary Yukawa potentials and they give an attractive 
contribution to the energy per particle. In the following we will refer this prescription as "without contact term". 
We will also consider a more realistic case where the pion and rho exchanges are folded with a two body correlation 
function G(r) As an example we take G(r) = jo(qc r) (qc = mj). This is equivalent to make a new change : 



q 2 + m 2 q 2 + m 2 q^ + m 2 

m 2 to 2 m 2 



q 2 + m 2 q 2 + m 2 q 2 -, + m 2 



(78) 



In conventionnal many-body language this is also equivalent to include a Landau- Migdal spin-isospin central interac- 
tion characterized by a g' parameter [l8|, [l9[ : 



Taking the pure VDM value, k p = 3.7, one obtains g' = 0.53, whereas the strong rho case, tt p = 6.6, leads to g' = 0.95. 
These values have to be compared with the most recent analysis which give g' ~ 0.6 [2(|. Since in this case it remains 
a contact term we will refer it as "with contact term" . Notice that the case "without contact term" can be recovered 
by taking the limit qc — ► 00; this corresponds to g' = 0.73 (k p = 3.7) and g' = 1.7 (n p = 6.6). 



A. Infinite symmetric matter 

We vary the ujNN coupling constant and the value of C in the (small) parameter region described above to 
reproduce the saturation density and the binding energy of nuclear matter : po = 0.16 fm -3 and E/A(po) = —15.96 
MeV. Another crucial parameter (although less well fixed) is the compressibility modulus K — 9p 2 d 2 e/dp 2 which 
should be around 250McV. As we will see, the present approach have difficulties to obtain such a low value for K and 
we choose to first reproduce the saturation density po and the binding energy E/A. 

Our results are given in Table 1 and the corresponding saturation curves are shown on fig. 2. For this discussion we 
keep the rho tensor coupling to its VDM value, namely k p = 3.7. We give our results with or without the dependence 
of kns on s and with or without the contact interaction. The first two lines of the table correspond to values of the 
parameters, (g^, C), giving a correct saturation point. In the last line of the table we give the asymmetry energy 
parameter 05, anticipating the next subsection. What come out from these calculations is a globally too large 
value of the compressibility modulus K even there is some numerical uncertainty in its determination and even if 
its extraction is not free of ambiguity in the context of relativistic theories [2l|. Removing the contact term {i.e., 
taking the Orsay prescription for the treatment of the pion and tensor rho exchanges) has the advantage to reduce 
the C parameter to a value closer to the lattice estimate, C ~ 1.25, but the incompressibillity K still increases by 
about 15 — 20 MeV. These contact terms in pion and rho exchanges (second term of the l.h.s of cq. [75)1 both give a 
repulsive contribution linear in density to the binding energy per nucleon: one finds 11.85 MeV for the rho contact 
piece and 0.6 MeV for the pion contact piece and the net attractive Fock term of the rho (—10.8 MeV) and the pion 
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Without contact 


With contact 






Kjvs(s) 




Hns(s) 




7.775 


7.678 


6.52 


6.42 


C 


1.33 


1.46 


1.49 


1.62 


p/po 


1.00 


1.00 


1.00 


1.00 


E/A(MeV) 


-15.96 


-15.93 


-15.97 


-15.92 


K(MeV) 


316 


297 


298 


281 


as (MeV) 


29.58 


29.45 


26.71 


26.64 



TABLE I: Values of the parameters and coordinates of the saturation point for different cases in symmetric matter. The last 
line is the asymmetry energy. The rho tensor coupling is k p = 3.7 (VDM value). The labels "without contact" and "with 
contact" interactions are explained in the text. The colums with kns(s) correspond to the case where the scalar susceptibility 
is density dependent and vanishes at full restoration. 




FIG. 2: Saturation curves for the cases including the contact terms with (solid) or without (dotted-dashed) the dependence of 
kjvs on s and cases without the contact terms with (dotted) or without (dashed) the dependence of kns on s. 

(—5 MeV) are thus reduced. Since these contact terms behave like an effective omega exchange, it is natural that the 
needed ujNN coupling constant is reduced when they are incorporated as apparent on Table 1 . The numbers relative 
to the Fock pion and rho exchange given just above correspond to the last column of Table 1 but they are actually 
very similar in the other cases. Incorporating the density dependence of the scalar susceptibility (kats(s) columns) 
has approximative^ the same effect as the absence of the contact terms in the sense that one has to increase the C 
parameter to compensate the lost repulsion. 

The main effect of the rearrangement terms is to reduce the effective nucleon mass by a few MeV but their influence 
on the saturation point is essentially negligible. This feature can be easily understood since these rearrangement terms 
enter the binding energy only through the dependence of the single particle wave functions on the Dirac effective mass 
M* = Mjv + £5 and in the non rclativistic limit this dependence just disappears. These rearrangement terms have 
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FIG. 3: Chemical potential with (solid line) or without (dashed line 
the energy per particle (dotted line) is also shown. 



inclusion of the rearrangement terms. For comparison, 



nevertheless an important influence on the chemical potential, which is defined by 

p=^=E*(k F ) + X°(k F ). 
dp 



(80) 



The Hugenholtz- Van-Hove (HVH) theorem |22j states that the chemical potential is equal to the binding energy per 
particle at the saturation point. On figure 3, we show the result of the calculation for the chemical potential p 
and for the binding energy per particle versus density, calculated again with the contact terms and with the density 
dependence of kns- We show the values of p with or without the rearrangement terms, the former being a few 
MeV smaller due the their effect on Ep. We see that in the presence of rearrangement terms, the theorem is exactly 
satisfied. Omitting the rearrangement terms yields to a significant violation of the HVH theorem. This failure has 
only a very moderate effect for the saturation curve of nuclear matter but would lead to major problems in finite 
nuclei, since the position of the Fermi energy would be displaced by about 5 MeV. 

As already mentionned, in the non relativistic limit the orbital wave functions are imposed: these are simply non 
relativistic plane waves. Consequently if the system is approximatively non relativistic, the results of the calculations 
should not depend very much on the choice of the wave functions. Said differently, the results obtained with another 
basis than the fully self-consistent Hartree-Fock (HF) basis would deviate from the HF results only by relativistic 
effects. For this reason we introduce the Hartree basis for which the Fock terms appearing in the self-energies (eq. 
171172173]) are removed. One advantage is the strong simplification of the calculations. The results for the saturation 
parameters are displayed on table 2. As expected, the results arc very similar but the results for C and for the 
asymmetry energy are slightly better in the Hartrcc-Fock basis. However the compressibility modulus decreases when 
the Hartree basis is used. 



B. Asymmetric matter 



The properties of asymmetric nuclear matter also constitute a crucial test for any effective theory. This is partic- 
ularly true for the study of nuclei far from stability and the equation of state of asymmetric matter is an important 
input in astrophysical studies and in particular for neutron stars properties. For that purpose let us introduce the 
asymmetry (isospin) parameter: (3 = (p n — p p )/p. 
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50 




FIG. 4: Asymmetry energy versus density. The calculation is done with the inclusion of the contact terms and when the scalar 
polarizability kns depends on s. 





Without contact 


With contact 




KNS 


KJVs(s) 


KNS 


Kjvs(s) 




7.728 


7.615 


6.45 


6.36 


c 


1.47 


1.62 


1.61 


1.735 


p/po 


1.00 


1.00 


0.99 


1.00 


E/A(MeV) 


-15.96 


-15.96 


-16.01 


-15.98 


K(MeV) 


301 


284 


286 


274 


a sym (MeV) 


26.15 


26.13 


24.10 


24.11 



TABLE II: The same as Table 1 but in the Hartree basis. 



a. Asymmetry energy. The first quantity to be considered is the asymmetry energy which can be defined according 
to : 



as{p) 



~d(E/A)~ 




\d{e/py 




dj3 


0=0 




(3=0 



(81) 



The contribution of the kinetic+Hartree term (eq. loT))) can be calculated analytically. In a pure Hartree approximation, 
omitting the extremely small contribution from the (5 meson exchange and the very tiny dependence of the scalar 
density ps on f3, one obtains the familiar result : 



(as) 



kin-\-Hartree 



F + 



M* 



2m 2 p 



In the actual calculation (last column of Table 1) the kinetic energy term gives (as)kin = 13.9 McV, and the Hartree 
rho exchange term gives (as) rho- Hartree— ^ ^07 McV. Hence The Hartree piece, (as)Hartree = 21 MeV, is not sufficient 
to reproduce the asymmetry energy. It follows that the Fock terms are absolutely necessary to get an higher value 
closer to the accepted values centered around 30 MeV. As for the binding energy there is a strong compensation 
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3.7 


5 


6.6 




6.42 


7.348 


8.62 


C 


1.62 


1.32 


0.93 


p/po 


1.00 


1.00 


1.00 


E/A(MeV) 


-15.92 


-16.05 


-15.94 


K(MeV) 


281 


313 


308 


a sym (MeV) 


26.64 


29.77 


35.87 



TABLE III: Evolution of the asymmetry energy with the value of the p-tensor coupling to the nucleon. The calculation is done 
with the inclusion of the contact terms and when the scalar polarizability k,ns depends on s. 



between the scalar and omega exchanges: (as) s +u = —0.75 MeV and the needed repulsion comes from the pion 
exchange {as) P ion = 2.75 MeV and mainly from the rho exchange (as)rho-Fock = 5.33 MeV. The important point 
is the specific contribution of the tensor piece of the rho exchange (as)rho-tensor = 8.03 — 4.12 ~ 4 MeV, the first 
number corresponding to the calculation within the Orsay prescription and the second (—4.12) coming from the 
contact term. The full result. 05 = 26.64 MeV, is nevertheless a little too low. However one has to keep in mind that 
the p-tensor piece goes like k 2 and that the p-tensor tensor coupling is not very well known phenomenologically. In 
the above discussion we took the pure VDM value. n p — 3.7, but the understanding of ir — N scattering data requires 
k p = 6.6 fl3j . Given this uncertainty we decided to vary this coupling between these two limiting values (See Tabic 
3). For K p — 5 we find as — 30 MeV. The new fit of the saturation point gives a value of the nucleon scalar response 
parameter, C = 1.32, very close to the lattice estimate but the weak point is that the compressibilty modulus increases 
beyond 300MeV. Independently of the details of the model calculation the above discussion strongly suggests that 
this p-tensor Fock term, which is absent in RMF theory, is a crucial ingredient to get a consistent description of 
bulk properties of symmetric and asymmetric nuclear matter. Its importance has also been stressed in the framework 
of density-dependent relativistic Hartree-Fock theory (DDRHF) where it is shown that it significantly improves the 
results for the single particle spectra in finite nuclei j23[. For completeness we also plot on fig. 4 the asymmetry 
energy against the density. We see a monotonic increase of as but with a slope decreasing at high density, which 
is very similar in shape to the result of a BHF calculation (24|. In any case this curve never exhibits the maximum 
found in the first version of in-medium chiral perturbation calculation (25| . 

b. Proton and neutron effective masses. The nucleon effective mass m* is a key property characterizing the 
propagation of the nucleons through the nuclear medium. In very exotic systems, the isosvector dependence of this 
effective mass, i.e., the mass splitting between the proton effective mass, m*, and the neutron effective mass, m* with 
increasing asymmetry f3 should play an important role. However, so far, no experimental data from finite nuclei has 
allowed a determination of the effective mass splitting but ab-initio Brueckner- Hartree-Fock calculation [24|, HfiJ have 
predicted m* > m* in neutron-rich matter. This result has been confirmed by relativistic calculations within the 
Dirac-BHF scheme [27| or within the DDRHF scheme (28|. Thus the sign of the splitting is rather solidly predicted 
although its amplitude is subject to a much greater uncertainty. This effective mass, which is a momentum dependent 
quantity, is defined from the density of states according to : 

to* k dk ^ 
where e is the single particle energy with the bare nucleon mass subtracted. In our approach it corresponds to : 



e = e£ - M N = J{M N + E s (k)) 2 + k 2 (1 + £y (k)) 2 + S (k) - M N 



With this definition, one finds : 



E* \I{M N + Ss (k)) 2 + k 2 (1 + (k)) 2 



(1 + 5V) 2 " (1 + Sy) 2 

Since there are many effective masses introduced in the literature a few remarks are necessary. 



(83) 



1. In the following we will discuss this effective mass taken at the Fermi momentum. This effective mass is called 
the Landau mass. 




2. This Landau effective mass should not be confused with the Dirac mass M* = Mn + Eg(k)) introduced 
above. The Dirac mass usually has an opposite behaviour to the Landau mass : in neutron-rich matter, 

M* <r M* 

lvl Dirac{n) ^ lvl Dirac(p) • 

3. It corresponds to the group mass in the terminology of ref. (29| and to the relativistic mass in the terminology 
of ref. [23I . Notice that the non relativistic mass M^ R introduced by these authors is not directly relevant for 
our discussion. 

Again we consider the case corresponding to the last column of Table 1 but the conclusions would be very similar in 
the other three cases. These effective masses arc plotted on fig. 4 which demonstrates that the behaviour characterized 
by m* > m* also occurs in our approach. It is also important to mention that the effect of the rearrangement terms 
is to decrease these two masses by about 7 MeV independently of the value of the neutron-richness parameter (3. 

c. Neutron matter and cold uniform matter in beta equilibrium. In this paragraph we determine the equation 
of state of cold uniform matter which might be relevant for studies of the interior of neutron stars. In this quite 
preliminary work we take a simplified point of view where this matter is mainly made of neutrons in beta equilibrium 
with protons and a gas of free electrons. Hence what is presented below should be considered only as indicative results 
in view of future more detailed work. The energy density of such a matter is therefore : 

£ ^nuclear matter ~i~ (^^) 

where e nuc i ear matter has been calculated before and e e - is 

^=1 E VAfT^cg| (85) 

fee <fcp e 



if we neglect the electron mass. To ensure the charge neutrality of the matter, we have to impose the equality between 
the densities of protons and electrons : p p = p e which implies pp = kp e . The grand potential and its density are 
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given by 



n = e - V Mi-^i 



c 



MS (Pp + pn) ~ Me (Pe - Pp) 



= £ - flpPp - H n p n - p e p e (86) 

wich implies ps — Pn = Pp + Pe- For a given total baryonic density, the proton and neutron densities and Fermi 
momenta can be obtained by solving the system of equations : 

kpe = Pf = e n {n F ) - e p {p F ) 

p = p p + p n . (87) 

For different values of p, we search a solution of this system around p p ~ O.f p and p n ~ 0.9 p. With those values of 
p p and p n , we can calculate the energy density versus baryonic density. We plot the resulting curve on fig. 6 for the 
case corresponding again to the last column of Table 1. We also determine the equation of state, namely pressure 
versus energy density. Two definitions of P can be used. The thermodynamical one is: 



P=-u; = -e + J2 



PiPi = —e + p n p 



where the charge neutrality condition p p = p e has been used. The other definition of the baryonic pressure is : 



2 5 (e/p) <9 £ 



P= 'i^f=^- e - (89) 

Since p is defined by ps = de/dp, the two definitions of P are thus equivalent. In the same conditions as above, 
we plot the pressure versus the energy on fig. 7. As all realistic calculations we find that the pure neutron matter 
is unbound and the energy per particle monotonically rises with the density. The equation of state (fig. 7) is rather 
similar to the one obtained with the QMC model [17| in the density range considered (p < 3po)- According to this 
latter work it is clear that beyond this density the contribution of the hyperons has to be included to make a realistic 
study of neutron star properties. 



VIII. CONCLUSION 



In this work we have investigated the implications for the nuclear matter problem of an effective chiral relativistic 
theory where the nuclear binding is ensured by a chiral invariant scalar field associated with the radial fluctuations 
of the chiral quark condensate. Nuclear matter stability can be obtained only if the response of the nucleon to the 
scalar field is properly incorporated, which phenomenologically produces three-body repulsive forces. The saturation 
mechanism depends on a balance between these repulsive forces and attractive three-body forces originating from 
chiral tadpole diagrams generated by the mexican hat effective potential. Once the various inputs are constrained 
by lattice QCD data and hadron phenomenology we obtain encouraging results with very minimal assumptions. One 
very positive point concerns the properties of asymmetric nuclear matter. We obtain the trend for the neutron-proton 
Landau effective mass difference to increase with increasing N — Z . Similarly the asymmetry energy, as, can be 
reproduced but a very important contribution comes from the tensor piece of the Fock rho exchange term. It is 
important to notice that the rho Hartree term is by far not sufficient if we keep the VDM value for the pNN coupling 
constant. One difficulty of the approach is however a too large compressibility modulus. This is one reason among 
other to motivate future works. The first line of investigation is related to the many-body treatment. Here we limited 
ourselves to the Hartre-Fock approach in the static approximation but higher order many-body effects, such as the 
inclusion of the correlation energy, have to be studied. For what concerns the pion loop correction this certainly 
necessitates to relax the static approximation. This inclusion of pion loops on top of the mean-field HF approach can 
be performed using in-medium chiral perturbation or many body approaches (RPA) in the line of our work of ref. 
Q but extended to the relativistic case. According to our non relativistic results [1] we expect that this correlation 
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P [fm 3 ] 

FIG. 6: Energy per particle versus baryonic density for cold uniform matter in beta equilibrium. The calculation is done in 
the scheme corresponding to the last column of Table 1. 

energy may decrease the compressibility modulus. Another line of research is related to the very crucial problem of 
the nucleon substructure. Here we have implicitly assumed that the nucleoli is made of three constituent quarks (or 
possibly a constituent quark and diquark) of the NJL type linked by strings. Although very attractive this picture 
has some problems. For instance it seems difficult to reproduce the value of the parameter kjvs characterizing the 
scalar response of the nucleon In such a picture the constituent quarks (or diquark) essentially move in the 
non pcrturbative vaccuum. However it may very well happen that the confinement mechanism generates a small 
region where the quarks are present with the highest probability. In such the region the chiral quark condensate (and 
consequently the quark mass) will drop due to the presence of the valence quark scalar density. In other words the 
confinement mechanism may generate a '"bag"' in which chiral symmetry is restored. We are thus back to the old 
dilemmna of the nucleon structure. Does the nucleon look like a chirally restored bag or does it look like a Y shaped 
string with constituent quarks at the ends? We believe that the scalar response of the nucleon or the chiral expansion 
parameter which are both relative to second derivatives of the nucleon mass are quantities which are very sensitive 
to the nucleon substructure. Certainly important efforts have to be done to obtain an estimate of these quantities for 
various modellings of the nucleon. 
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FIG. 7: Equation of state of cold uniform matter in beta equilibrium : pression versus energy density. The calculation is done 
in the scheme corresponding to the last column of Table 1. 

APPENDIX 



APPENDIX A: EXPLICIT FORM OF THE FOCK TERM ENERGY 

The various contributions to the Fock term energy, completing cq. 1441 arc listed below: 



E 



Fock 



j J dvdr' Tr(S p (v' - r) S p (r - r') + S n (r' - r) S n (r - r')) 
<£(r)<£(r')A,(r-r') 



(Al) 



E 



Fock 



J drdr'Tr(S p (r'-r)^S p (r-r') 7 v + 5„(r' - r) 7 /1 5„(r - ^7")) 
D^ v {r - r') 



(A2) 
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E fL = -f J dvdv'Tr(S p {v' -v)YS p (r~v'h 1 ' + S n (r' - r)^ S n (r - v')j v 

+ 2S p (r'-r) 7 ^(r-r')7" + 2S n (r' - r)f5 p (r - r') 7 v )) 

" 2 f (jM^) J drdr ' Tr (S P (r'-r)^S p (r-r')^ 
lS„(r'-r)a«S„(r-r'y 

+ 2S p (r' - v)a^S n (r - r')Y + 2S„(r' - v)a^ S p (r - v')j v )) 



2 / \ 2 r 

/ drdr' Tr(S p {r' - r)a^S p (r - r')^ v 



2 \2M 



N 



+ S„(r^rK'S„(r-r')/' 

+ 2 S p (r' - r)<j^S„ (r - r')a v i + 2 S„ (r' - r)a^S p {v - r')a uj ) 
did'jDp^ir-r') (A3) 

E fL = f j ^ ^ ^ ( Sp( - T ' ~ ^ 5p(r " + 5 " (r ' " ^ ^ ~ 

+ 2 5 p (r'-r)5„ (r-r') + 2 S n (r' - r) S p (r - r'))) D s (r - r') (A4) 



41* = ^(f£) / drdr'Tr^ p (r'-r) 7 V^(r-r')7V 

+ 5„(r'-r) 7 5 7 'S n (r--r')7V 
+ 2S p (r'-r)7V^(r-r')7V 

+,2 S n (r' - r) 7 5 7 4 S* p (r - r') 7 V )\ O^D^r - r') 



(A5) 
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APPENDIX B: EXPLICIT FORM OF THE ENERGY DENSITY IN INFINITE NUCLEAR MATTER 



The various pieces of the energy density, completing cq. (|70p . are listed below. For simplicity the occupation 
numbers, N^N^' , arc omitted. 



W _ 9 S 



Fock 



e 



c 



ir VT 

c 



Sp 1 



f f d 3 k d 3 k' 1 ^ M* M'* k* k 



El XVI XVI 14. 14. \ . . 

1 + ( B1 ) 



2 7 (2tt) 3 (2tt) 3 (k - k') 2 + m* 2 ^ ^ E*E'* E* E'* J 
Focfc " 2 7 (2tt) 3 (2tt) 3 (k - k') 2 + ml ^\ E* E'* E* E'* J 



M _ 9s 

"Fock 



2 f d 3 k d 3 k' 1 , , / M* M'* k* k' 



r d 3 k d 3 k' 1 ^ / 

J (2^) 3 (2tt) 3 (k-k'^+m 2 ^! 1 



2 J (2tt) 3 (2tt) 3 (k - k') 2 + m\ ^ ^ £>£T* £> 

/ ^ 3k d 3 ^ i A , W Ml _ ^ ™ 
+ 55 y (2^) 3 (2tt) 3 (k - k') 2 + mi \ + £5 k;* e; ' k* ; 1 ] 

(p v ) 9 2 P f d 3 k d 3 k' 1 ( n t M* M'* „k* k'*' 



E XVI XVI K. 

2 /■ d 3 k d 3 k' 1 / Af*M^* k; k^*\ 

- 29 "J (2ir) 3 (2tt) 3 (k - k') 2 + m 2 I 2 " ~ ^ ' K* ) { ' 



CP'"") = 6 1 /, /• ^ ^ 1 ^ / k*M*\ 

Focfc 2 yp 2M N J (2tt) 3 (2tt) 3 (k - k') 2 + m 2 ^ ^ ; E* E'* j n 

1 / p /■ d 3 k d 3 k' 1 



r p 2M N J 



(2tt) 3 (2tt) 3 (k-k') 2 +m 2 



X I fk-k't • ^L^L. + Ck- k') • ^k-^L | (B5) 
1 / f p \ 2 f d 3 k d 3 k' 1 2 M*M'* 



2\2M N J J (2tt) 3 (2tt) 3 (k - k') 2 + m 2 ^ V ' E*E'* 

.( k - k') 2 — ■ ^ +4(k - k') • — (k- k') • ^ 

/ /„ \ 2 r d 3 k d 3 k' i a ,,. 2 , Q .. ^. 2 m;m'* 



V2M N ; y (2^) 3 (2^) 3 (k-k') 2 + m 2 1 ( j ; 



p n 



-(k - k') 2 § ■ ^ + 4(k - k') • §(k - k') • 4) (B6) 

e W _ IfM 2 /" d3k d " k ' 1 Vf( k kQ 2 I (1: U)i ArM '* 

Fock ~ 2 \mj J (2tt) 3 (2^) 3 (k - k') 2 + ml ^ [} ' \ _E* £7'* 

_( k - k') 2 — • iV + 2(k - k') • — (k - k') • 

v ' E* E * v ' E* v 7 £ * 



J (27r) 3 (27r) 3 (k-k') 2 + m 2 ^ k j + (k kj 



m^y y (2tt) 3 (2tt) 3 (k - k') 2 + m 2 V ^ E'* 

< k k ') 2 § ■ ft + 2(k - k') - |§ (k - k') • ft) (B7) 



27 



B.D. Serot, J.D. Walecka, Adv. Nucl. Phys. 16(1986) 1; Int. J. Mod. Phys. E16 (1997) 15. 

P. Finelli, N. Kaiser, D. Vretenar and W. Weise, Nucl. Phys. A770 (2006) 1. 

W. Bentz and A.W. Thomas, Nucl. Phys. A696 (2001) 138. 

G. Chanfray, M. Ericson and P.A.M. Guichon, Phys. Rev C63 (2001)055202. 

M. Ericson and G. Chanfray, EPJA 34 (2007) 215. 

A.K. Kerman and L.D. Miller in "Second High Energy Heavy Ion Summer Study", LBL-3675, 1974. 

G. Chanfray and M. Ericson, EPJA 25 (2005) 151. 

G. Chanfray and M. Ericson, Phys. Rev C75 (2007) 015206. 

T. Hatsuda, T. Kunihiro and H. Shimizu, Phys. Rev. Lett. 82 (1999) 2840. 

A. W. Thomas, P. A. M. Guichon, D. B. Leinweber and R. D. Young Progr. Theor. Phys. Suppl. 156 (2004) 124; 
|nucl-th/04fT014l 

P.A.M. Guichon, Phys. Lett. B200 (1988) 235. 

F. Bissey et al, Nucl. Phys. B (Proc. Suppl.) 141 (2005) 22, |hep-lat/0501004| 



G. Hohler and E. Pietarinen, Nucl. Phys. B95 (1975) 210. 
P. A. M Guichon, H. H. Matevosyan, N. Sandulescu and A. W. Thomas, Nucl.Phys. A772 (2006) 1. 
P. Bernados et al, Phys. Rev. C48 (1993) 2665. 
A. Bouissy et al, Phys. Rev. C36 (1987) 380. 

J. Rikovska-Stone, P. A. M Guichon, H. H. Matevosyan and A. W. Thomas, Nucl. Phys. A792 (2007) 341. 
W. Weise, Nucl. Phys. A278 (1977) 402. 

G. Chanfray, J. Delorme and M. Ericson, Phys. Rev. C31 (1985) 1582. 

M. Ichimura, H. Sakai and T. Wakasa, Prog. Part. Nucl. Phys., 56 (2006) 446. 

G. Colo and Nguyen Van Giai, Nucl. Phys. A731 (2004) 15. 

N. M. Hugenholtz and L. Van Hove, Physica (Amsterdam) 24 (1958) 363. 

W-H. Long, H. Sugawa, Nguyen Van Giai and J. Meng, Phys. Rev. C76 (2007) 034314. 

W. Zuo, J. Bombaci and U. Lombardo, Phys. Rev. C 60, 024605 (1999). 

N. Kaiser, S. Fritsch and W. Weise, Nucl. Phys. A697 (2002) 82. 

J. Bombaci and U. Lombardo, Phys. Rev. C 44 (1991)1892. 

F. Sammarruca, W. Barredo and P. Krastev, Phys. Rev. C71 (2005) 064306. 

W-H. Long, Nguyen Van Giai and J. Meng, Phys. Lett. B640 (2006) 150. 

M. Jaminon and C. Mahaux, Phys. Rev C40 (1989) 354. 



o 
o 



o 
o 

(N 



o 
in 



o 
o 



o 
in 



o 
o 



O 
O 

in 



o 
o 



o 
o 

m 



o 
o 

(N 



o 
o 



